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Abstract 

The present study numerically demonstrates how 
small spanwise variations in velocity upstream of 
a body can cause relatively large increases in the 
spanwise- averaged heat transfer to the leading edge. 
Vorticity introduced by spanwise variations, first de- 
cays as it drifts downstream, then amplifies in the 
stagnation region as a result of vortex stretching. 
This amplification can cause a periodic array of three- 
dimensional structures, similar to horseshoe vortices, 
to form. The numerical results indicate that, for 
a given wavelength, there is an amplitude thresh- 
old below which a structure does not form. A one- 
dimensional analysis, to predict the decay of vorticity 
in the absence of the body, in conjunction with the 
full numerical results indicates that the threshold is 
more accurately stated as a minimum level of vortic- 
ity required in the leading edge region for a structure 
to form. It is possible, using the one-dimensional anal- 
ysis, to compute an optimum wavelength in terms of 
the maximum vorticity reaching the leading edge re- 
gion for given amplitude. A discussion is presented 
which relates experimentally observed trends to the 
trends of the present phenomena. 

Nomenclature 

A peak to peak velocity difference 

normalized by the average velocity 
d diameter of wires in screen 

Fr Frossling number (Fr = Nu/y/Re) 

k thermal conductivity 

M Mach number 

Nu Nusselt number 

p pressure 

Qw wall heat transfer 

R leading edge radius 

Rcj; Reynolds number (i.e., Re 3; = poUx/po ) 

T , temperature 
u streamwise velocity 

V normal velocity 

* Member AIAA 


w 

spanwise velocity 

Ux 

velocity difference (ie. u\ = u — U ) 

^mtn 

minimum streamwise velocity in 


X — z symmetry plane 

u 

spanwise- averaged freestream velocity 

X 

streamwise direction 

Ax 

distance from xq 

^wire 

denotes position of screen 

y 

normal direction 

z 

spanwise direction 

Greek Symbols 

A 

wavelength of spanwise variation 

T 

ratio of specific heats 

A* 

viscosity of fluid 

n 

vorticity oriented in the y-direction 

fi 

normalized vorticity predicted by the 


analysis (ie. Q == Q/{7tU /R)) 

p 

density of fluid 

a 

solidity of array of wires 

e 

distance along surface normalized by the 


leading edge radius 

Subscript 

.s and Superscripts 

Oe 

conditions at the exit 

Of 

represents the mean value in the 


boundary layer (ie. film) 

OC.B 

denotes leading edge 

()rm'n 

represents the minimum value on the 


stagnation line 

Omaj" 

represents the maximum value on the 


stagnation line 

Oopt 

wavelength for which analysis predicts 


maximum peak vorticity at leading edge 

Ot 

total conditions 

OtH 

threshold point at which structure forms 

Ou- 

value at wall 

()ii;*re 

denotes position of array of wires 

Oo 

denotes upstream boundary 

0 

span wise- aver aged quantity 
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Figure 1: Geometry of Problem 


I Introduction 

The prediction of heat transfer to the leading 
edge of gas turbine blades is very important to the 
design process. Typically, the leading edge is a lo- 
cation of high heat transfer. It has been known for 
some time tliat relatively small increases in freestream 
turbulence above a nominally laminar flow can cause 
large increases in the heat transfer at a forward stag- 
nation region [1]. However, the mechanism by which 
such an increase occurs is not well understood. It has 
also been observed experimentally that span wise varia- 
tions in the heat transfer occur, with wavelengths that 
vary inversely as the square root of the Reynolds num- 
ber, in flows with very small spanwise variations in the 
freestream velocity. This phenomena has been studied 
experimentally as it influences heat transfer [2, 3] and 
mass transfer [4, 5, 6]. The experiments indicate that 
spanwise velocity variations on the order of 1% can 
cause the transfer coefTicients to increase on the order 
of 20%. 

Asymptotic analysis for small perturbations in the 
spanwise direction [7, 8, 9, 10, 11] indicate that distur- 
bances with a wavelength greater than some neutral 
wavelength can be amplified in the boundary layer. 
However, prediction of the neutral wavelength and the 
most unstable wavelength differs from author to au- 
thor. In fact, Hammerlin [12] found instability for a 
continuous spectrum of wave numbers, while Wilson 
and Gladwell [13] found the flow was always stable 
to infinitesimal disturbances. More recently, Thomp- 
son [14] found that stagnation flow exhibits instability 
if unsteadiness is introduced into the basic state of the 
flow. His analysis suggests that there is an amplitude 
threshold on the unsteadiness for instability to occur. 

la a revdew of this type of flow, Morkovin [15] 
discusses the limits of the available analysis tools as 
well as the uncertainties in experiments. It is unclear 
whether the major effects are a result of the random 


unsteady turbulence or a result of large, essentially 
steady, three-dimensional structures. These structures 
could be the result of small spanwise perturbations in 
the freestream. Morkovin [15] also raises the ques- 
tion of forced or free response: “The general impres- 
sions seem to be that some forcing is present and tliat 
somehow the characteristics of the boundary layer on 
the cylinder enter the quantitative determination of 
the geometry of the horseshoe vortices.” Based on a 
combination of flow visualization, turbulence measure- 
ments, and spanwise heat transfer measurements, Van 
Fossen and Simoneau [2] arrived at the same general 
impressions. 

Numerical studies attempting to model the tur- 
bulence [16, 17, 18] have also been undertaken with a 
certain degree of success. However, if the experimental 
results are dominated by a steady spatial variation or 
even a nonrandom variation in the third dimension, a 
two-dimensional calculation with a turbulence model 
will not properly predict the flow features. 

The present work will concentrate on the effect of 
a steady spanwise velocity variation on the heat trans- 
fer at the leading edge. Numerical solutions of the full 
three dimensional Navicr-Stokes equations are used to 
investigate the requirements for the formation of a 
horseshoc-like vortex, and to show the effect of these 
vortices on the heat transfer coefficient at the leading 
edge. In addition, a parallel flow analysis which pre- 
dicts the decay of v’orticity in the absence of the body is 
employed. This analysis proves to be useful in predict- 
ing the threshold level of vorticity for the formation of 
a vortex structure. The parallel flow analysis also pre- 
dicts the optimum wavelength in terms of supplying 
the maximum vorticity to the leading edge region. 

II statement of Problem 

Consider the problem of a flat plate with a cylin- 
drical leading edge of radius R, Figure 1. The body 
is assumed to be infinite and uniform in the spanwi.se 
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Figure 2: 37x51x22 grid for A = 0.4/? 


direction, and semi-infinite in the streamwise direc- 
tion. Upstream of the body, the flow has a sinusoidal 
variation of velocity in the spanwise direction with a 
wavelength A, which introduces vorticity oriented per- 
pendicular to the X -- z plane. The heat flux on the 
body is fixed at a uniform value. The peak to peak 
velocity variation as a percentage of the average up- 
stream velocity is defined as Aq and the wavelength is 
denoted by A. The velocity distribution is imposed at 
a position upstream of the leading edge, denoted by 
xq. This position is very important to note since the 
vorticity first decays as it approaches the body, then 
is amplified in the region near the leading edge. By 
ignoring the presence of the body, the rate of decay 
of the vorticity can be estimated from a parallel flow 
analysis, as will be shown below. 

Ill Solution Method 

The solution was computed using the PARC3D 
code, which uses the Beam and Warming approxi- 
mate factorization algorithm, a time marching im- 
plicit scheme that solves the Navier-Stokes equations 
in strong conservation form. It uses second order 
accurate central- differencing for spatial derivatives. 
The code was originally developed as AIR3D by Pul- 
liam and Steger [19]; Pulliam [20] later added the 
Jameson [21] artificial dissipation and called the code 
ARC3D. Cooper [22] adapted the code for internal 
flows in propulsion applications and named the code 
PARC3D. 

For the present application, the equations are 
solved for the steady laminar flow solution in the 
leading edge region, with the algebraic Baldwin- 
Lomax [23] turbulence model turned on just down- 


stream of the juncture between the cylinder and flat 
plate. During the initial part of the calculation, the 
turbulence is updated at each iteration; then it is 
frozen and only recalculated each time the calculation 
is restarted, roughly every 2500 iterations. This freez- 
ing of the turbulence conserves computer time, and 
enhances the rate of convergence. 

The turbulence model is used here simply to en- 
hance convergence by limiting any reverse flow wliich 
may occur in this region. The primary focus is the 
leading edge region, which should not be greatly af- 
fected by the turbulence downstream. Because of the 
central-differencing, the code has Jameson [21] type, 
fourth order, artificial dissipation for stability. The ar- 
tificial dissipation is added both explicitly and implic- 
itly; the coefficient varies both spatially and direction- 
ally to minimize the dissipation [24, 25]. Although a 
dissipation model which varied the coefficient spatially 
but not directionally was too dissipative, it gave very 
good convergence rates. Because of this, the local- 
directional model WcLS used both upstream of and on 
the cylindrical leading edge, and the local model was 
used both on the flat plate portion and downstream 
of the leading edge. In a problem of this type, it is 
especially important to use a method where the de- 
cay of frees tream vorticity from numerical dissipation 
is small compared to the decay from the physics of the 
flow. Included in the results is a demonstration that 
the decay of freestream vorticity is properly calculated. 

A C-grid, shown in Figure 2, was used; half of the 
grid was eliminated by imposing symmetry in the x — z 
plane. The grid extended 10 radii up and downstream 
and ranged from z = A/2 to z = A. 
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The PARC3D code requires about 30 words of 
storage per grid point and 22 x 10*"^ seconds of CPU 
time per iteration per grid point on the NAS Cray 
Y/MP. Each case took approximately 5000 iterations, 
for typical run times on the order of 5000 CPU sec- 
onds. For each value of A, the problem was first solved 
for some value of Aq, then Aq was changed, and the 
previous solution was used as an initial guess. This re- 
sulted in fewer iterations required for subsequent runs. 

Boundary Conditions 

A symmetry condition was imposed at z = A/2 
and z = A. Symmetry was imposed by using second 
order accurate one-sided differences. 

On the body, no-slip, constant heat flux conditions 
were enforced, along with a zero normal pressure gra- 
dient. At the exit, static pressure was imposed while 
velocity components and total energy were linearly ex- 
trapolated. The upstream boundary condition is de- 
termined by specifying the total temperature and total 
pressure, and requiring that the flow be aligned with 
the x-direction. Then at each iteration, the upstream 
running Riemann invariant is extrapolated to predict 
the velocity and density. The original code used ze- 
roeth order extrapolation, evaluating the invariant one 
grid point in from the boundary and imposing it at the 
boundary. Using this technique resulted in a ”kink” 
in the decay of vorticity at the upstream boundary; 
that is, the vorticity showed very little decay between 
the first two grid points and then decayed normally. 
This problem was greatly reduced by using a linear 
extrapolation for the Riemann invariant. 

To produce a sinusoidal velocity variation at the 
upstream boundary, the total pressure was specified 

by 


+ ( 1 ) 

where M' is approximately the peak Mach number at 
the upstream boundary and Aq will be roughly the 
peak-to-peak velocity difference as a percentage of the 
inlet velocity. This boundary condition will produce a 
Mach number and amplitude slightly different from A/' 
and Aq. The results are reported based on the actual 
solution, not on the idealized parameters. For all of the 
present cases ^ 0.2, resulting in freestream Mach 
numbers of roughly 0.17. The total temperature was 
specified by 

T, =To{l + ^[A/'(l - ^(1 + co5(^))]='} (2) 

This total temperature profile was imposed simply as 
an aid to assess the numerical dissipation of the code. 

This profile produces a uniform entropy profile at 
the upstream boundary. By running the code without 


the viscous terms, it was possible to see the effects of 
the numerical dissipation in terms of the production of 
entropy. The inviscid calculations were used to indi- 
cate the number of grid points required for an accurate 
calculation of the decay of vorticity in the freestream. 
The total temperature profile was left as stated above 
for all the runs, since its effect is negligible at the low 
Mach numbers. 

IV Decay of vorticity in the 
freestream 

To approximate the rate of decay of vorticity be- 
tween the upstream boundary and the body, consider 
the flow of an incompressible fluid in the x ^ z plane. 
Flow is assumed to be in the x-direction only, with 
a periodic variation in the z-direction. In addition, 
assume that the pressure gradient in the x-direction 
is zero and that the initial periodic variation is small 
compared to the average velocity. Then, let the veloc- 
ity difference be 


Ui(x,z) = U - u{x,z) (3) 


where U is the average velocity. The appropriate equa- 
tion governing the velocity difference, ui, is 


du\ ^ u d^ui 
dx U dz^ 


(4) 


assuming Ui/tJ << 1, where the boundary layer as- 
sumption has been made and terms quadratic in ui 
and Wi have been neglected; Wi is the cross stream 
velocity, which is being neglected. The boundary con- 
ditions are 


at X = xo 
periodic 


AoU 27TZ 

ui(xo, W 

wi(a;,0) = «i(a:, A) (6) 


where Aq is an arbitrary constant, but must be small 
compared to one. The solution to this equation with 
these boundary conditions is 


ui(x,z) 


AqU 

2 


cos(— — )exp(— 

A 


47T^ (x - Xq) 
Rex A 


) (7) 


where Re\ = UXfi' . Of particular interest is the 
vorticity, which is oriented out of the x — z plane, and 
can be written as 


n(x,z)^^ 


( 8 ) 




Notice that the level of vorticity is everywhere directly 
proportional to the constant Aq, but the rate of decay 
is strongly dependent on the wavelength. Increasing A 
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Figure 3: Convergence History for A = 0.4/2, A = 0.04, Rcr = 6500 


decreases the initial value of vorticity at z = xq, but 
it also decrecises the rate of decay so that at some dis- 
tance downstream, the level of vorticity will be greater 
than that for a case with a smaller A and the same Aq. 
The wavelength which results in the maximum peak 
vorticity at some distance downstream, to be referred 
to as the optimum wavelength X^pti can be written as 

where Ax - x - Xq and /2eAx = lJ{x ~ xo)/u . Note 
that the optimum wavelength is directly proportional 
to y/ Ax for a given velocity, and inversely proportional 
to y/Re^j; for a given Ax. 

The main result of this analysis is the finding that 
an upstream variation with the same amplitude but a 
longer wavelength will not necessarily create a smaller 
disturbance at the leading edge of the body. Further 
discussion of this point will appear below as it pertains 
directly to the numerical results. 

V Results 

For the present study, the effect of varying the 
wavelength A and the amplitude Aq of the upstream 
velocity profile are considered. The Reynolds num- 
ber based on the leading edge radius is 6500, and the 
fixed wall heat transfer is set at = ikoTo/R for 
all cases, where ko and Tq are the freest ream ther- 
mal conductivity and temperature, respectively. This 
level of heat transfer results in a temperature rise of 
roughly 10% above the freestream absolute temper- 
ature. Calculations were run for three values of A, 
A = 0.2/2, 0.4/2, and 0.8/2. The parameter Aq varied 


through the range 0.0 to 0.066. The grid used for these 
cases was 37x51x22, in the streamwise, normal, and 
spanwise directions, respectively. The case of /Iq = 0 
was calculated using the same x — y grid as the three 
dimensional cases, but only five grid points in the z~ 
direction. The case of A = 0.2/2 and Aq = 0.066 was 
also calculated, using a grid of 76x101x42 to check 
grid dependency. For the cases with A = 0.2/2, uni- 
form grid spacing in the z-direction was used. After 
running the A = 0.4/2 Ccises, the decision was made 
to use geometric stretching in the spanwise direction 
to cluster points in the region where the vortex forms. 
The A = 0.4/2 cases were then rerun. Comparison of 
these two sets of runs is used as an additional check on 
grid dependency. The A = 0.8/2 cases were done with 
geometric stretching in the ^r-direction. 

Convergence Criterion 

For each case the code was run until the resid- 
ual dropped below 5 x 10“^®, corresponding to a drop 
of roughly three orders of magnitude from the initial 
level. The residual is defined as the sum of the L2 
norms for each of the five equations divided by the 
number of points. Figure 3 shows the history for the 
case of A = 0,4/2 and Aq = 0.04. This calculation was 
initialized using the converged solution for Aq = 0.026 
so the residual starts out relatively low then, and drops 
2.5 orders of magnitude to meet the convergence cri- 
terion, Figure 3a. The spikes in this figure result from 
the update of the turbulence each time the calcula- 
tion is restarted. Recall that turbulence is only cal- 
culated on the flat plate portion of the body and is 
frozen after the initial drop in the residual. Notice 
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Streamwise Distance, x/R 

Figare 4: Decay of Peak Vorticity, Aq = 0.04, Rcji ~ 
6500 

til at the span wise- averaged Frossling number at the 
leading edge, Figure 3b, has settled long before the 
calculation is terminated. The FYossling number is de- 
fined as the Nusselt number divided by the square root 
of the Reynolds number. This figure is representative 
of the other cases run. 

Decay of Vorticity in the Freestream 

Figure 4 demonstrates the decay of vorticity in 
the streamwise direction and the amplification near 
the leading edge. Shown in this plot is the y com- 
ponent of vorticity as a function of x along the line 
y 0, ? = 3A/4 . For the same value of = 0.04, 
the A = 0.2/?,0.4fl, and 0.87? results are shown. Also 
shown in this plot are the analytical results from equa- 
tion 9 with Ao defined so that the vorticity corresponds 
to the numerical result at x = Recall that the 
analysis does not take into account the presence of the 
body, and thus does not predict the amplification near 
the leading edge. This plot demonstrates the impor- 
tance of the rate of decay in determining the supply of 
vorticity to the leading edge region. Notice that the 
A = 0.27? case has the largest initial vorticity and yet 
does not result in the strongest vortex because of the 
high rate of decay. This plot also indicates that there 
is an optimum wavelength in terms of supplying the 
most vorticity to the leading edge for a given Aq. 

The optimum wavelength for the present Reynolds 
number and upstream location from equation 10 is 

A,,^, == 0.337? (11) 

which is in agreement with the observation, from the 
numerical computation, that A = 0.47? results in a 
stronger structure than the A = 0.27? or the A = 0.87? 



Figure 5: Decay of Normalized Peak Vorticity, Rcji — 
6500 

cases. Notice from equation 10 that t his optimum 
wavelength varies inversely with If Rcfi were 

varied while holding 7?/ Ax constant, then the previous 
statement would imply that the optimum wavelength 
varies as the inverse of y/Ren. It is very important 
to notice, however, that the optimum wavelength does 
not vary witli 7?, (i.e., if Reji is varied by changing 7?, 
then the optimum wavelength would not vary). Var- 
ious experiments [26, 27, 28, 29] with levels of turbu- 
lence on the order of 0.02 and less have observed span- 
wise variations with wavelengths that vary inversely as 
the square root of 7?e/^. This observation is often used 
to support the argument that the disturbance must be 
caused by the instability of the leading edge boundary 
layer. However, if the upstream variation is made up of 
many wavelengths, then the wavelength of the result- 
ing disturbance could be seen to vary inversely with 
y/Rcfi, when the freestream velocity is varied, because 
the optimum wavelength predicted by the parallel flow 
analysis varies 2 is y/Reji. It is not possible at this point 
to reach a conclusion, but it would appear that further 
investigation of the question of supply versus stability 
is warranted. 

Figure 5 demonstrates further the agreement be- 
tween the calculation and the analysis in the region 
of decaying vorticity. In this figure there are fifteen 
separate solid lines corresponding to ^ of the cases 
calculated; six cases with A = 0.27?, six cases with 
A = 0.47?, and three cases with A = 0.87?. By plot- 
ting the vorticity normalized by its upstream value, the 
common rates of decay for each wavelength at several 
values of Aq can be seen. Also shown in Figure 5 is the 
analytical result for each wavelength. The agreement 
is quite good considering the simplicity of the analysis. 
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Figure 6; Extrapolation to Threshold Level 


Once again, recall that the analysis does not take into 
account the presence of the body, and thus does not 
predict the amplification near the leading edge. 



It was observed numerically that a minimum am- 
plitude at each wavelength is required for a vortex to 
form. A vortex is considered to have formed when 
reverse flow occurs in the a: -* z symmetry plane. In 
Figure 6a, the minimum streamwise velocity Umtn is 
shown for three levels of Aq for each wavelength. The 
levels of Ao shown correspond to the first three values 
of Aq, which resulted in reverse flow for each wave- 
length. The predicted value of Aq for = 0 is then 
extrapolated by fitting a parabola thru the available 
data points. Figure 6a demonstrates that the mini- 
mum value for Aq which generates a vortex structure 
differs by roughly ±20% of the average for these three 
wavelengths. The peak level of vorticity supplied to 
the leading edge, as predicted by the parallel flow anal- 
ysis, is evaluated for each case from 


= ( 12 ) 

wU/R ^ ’ 

. R , 47t^ 2 (/? - a;o). 

= a) Se > 


a structure to form could be estimated from 


Aio,.. = (M) 


where ^L.E.,th is a constant and is deduced from previ- 
ous cases, and (/? — xq) is the distance from the leading 
edge to the upstream boundary. It is very important 
to recall that the present analysis is only valid in re- 
gions where Aq << 1; thus, [R — xq) could not be 
interpreted as the distance to a screen upstream, al- 
though changes in (/?— xq) would correspond directly 
to changes in the screen position. Notice that as the 
Reynolds number approaches infinity, the level of Aq 
required at the upstream boundary decreases to a con- 
stant. Averaging the three extrapolated values gives 
^L.E.,th = 0.037. Using this value in equation 14 pre- 
dicts the threshold level of Aq for each wavelength to 
within ±5% of the extrapolated values in Figure 6a. 

It is interesting to note that if the vorticity up- 
stream of the body occurs at the optimum wavelength 
as defined by equation 10, then the threshold level for 
Aq could be written as 

, , l{R-xo) 87t2 ^ 1 ^ 

Ao.tH\opi-d ^ 


where (R — xo)/i? = 9 for all of the present cases. 
This quantity is then extrapolated quadratically to the 
point of tirnm = 0. The predicted level of vorticity 
to the« leading edge differs by less than ±5% for the 
three wavelengths considered in Figure 6b. From this 
result it seems possible that, for a given wavelength 
and Reynolds number, the minimum Aq required for 


This result implies that if vorticity is always sup- 
plied at the optimum wavelength, the threshold Aq 
w ould vary in versely with y/Reji and directly with 
y/{R ^ ^o)/ R‘ 

Experimentally, the threshold level has been stud- 
ied by Hodson and Nagib [3] and by Kottke [4]. In each 
case, a row of wires was placed upstream of a circular 
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leading edge to produce a small spanwise variation in 
the upstream velocity. Kottke found that for Reynolds 
numbers such that the wire wakes were steady, the 
upstream position of the wires for impending vortex 
formation depended linearly on the Reynolds number 
based on the wire diameter and did not depend on the 
Reynolds number based on the leading edge diameter. 

It is possible to explain this trend in terms of the 
supply of vorticity reaching the leading edge. If the 
position xo is considered to be a fixed distance down- 
stream of the wires located at, say, then 

Ft Xq — It, Xu;i>e (^0 ^wire^ 


so changes in {R — xq) would correspond directly to 
changes in the distance from the wires to the leading 
edge. Requiring that the peak vorticity reaching the 
leading edge be a constant, equation 9 indicates that 
(7Z “ Xo) must vary linearly with Rex. This is equiv- 
alent to Kottke’s observation, since the ratio of wire 
spacing to wire diameter is not changed in that experi- 
ment. This conclusion also assumes that the dominant 
wavelength and Aq do not vary with Reji. Hodson and 
Nagib [3] found a correlation between y/R — x^ir^jd 
and (1 — (r)y/Reji, where a is the solidity of the screen 
and d is the wire diameter. They found a correlation 
for a single wire as well as for an array of wires with 
various wire diameters. The trend for each case agrees 
with the arguments of the present hypothesis. Ignoring 
the effects of solidity, tliis correlation is equivalent to 
a linear relation between {R — x^j^ire) and Rcr. Since 
Rcji was varied by changing the freest ream velocity, 
Rex will differ from Re^ by only a multiplicative con- 
stant. In fact, if a set of data does exhibit this linear 
relation except for some scatter from error, then tak- 
ing the square root of each quantity would tend to 
enhance the correlation. 


Heat Transfer Results 


Figure 7 shows the spanwise- aver aged Frossling 
number Fr as a function of 0 for ^ 0.04 and 
\/R — 0.2, 0.4, and 0.8, Also shown are the analytical 
solution of Frossling [30] and the PARC3D calculation 
for Aq — 0. The Frossling number is defined as 


Fr = 


Nu 


(17) 

(18) 


so 


= ^ [El(Zi!L_i] 

koTf o y/2ReR kj y /iq \Ti^a J 


(19) 


where is the Nusselt number, q^, is the wall heat 
transfer, Tu, is the wall temperature, and Tt^u is the 



Figure 7: Frossling Number vs. Streamwise Distance, 
Aq ~ 0.04, Rcfi ~ 6500 

freestream total temperature. The thermal conductiv- 
ity k and viscosity /i are evaluated at the film tem- 
perature to minimize the effect of variable properties. 
Their effect for the present work is not large, but 
it should be noted that evaluating properties at the 
freestream temperature can result in Frossling num- 
bers one to two percent higher than those presented. 
For the three-dimensional problem, the film tempera- 
ture is defined by 

Tj = 0.25(27f 0 + T^,max + (20) 

Notice that the Frossling analysis and the two- 
dimensional calculation (i.e., Aq = 0) show' good agree- 
ment for 0 < 30®. For 6 > 30® the Frossling analysis 
decreases faster than the two-diraensLonal calculation. 
This decrease is indicative of a comparison between 
a constant wall temperature condition, which is what 
the Frossling analysis uses, and a constant w'all heat 
flux condition [31]. Also notice in Figure 7, that for the 
same value of Aq, A//? = 0.4 causes the highest Fr, 
which is 18% higher at 0 = 0 than the two-dimensional 
result from PARC3D. In fact, the spanwise-averaged 
Frossling number was higher than the two-dimensional 
result for every case considered. Increases in the trans- 
fer coefficient on the order of 15 -- 20% have been seen 
experimentally [2, 3, 6] in flows with upstream span- 
wise variation. Direct comparison with the present re- 
sults would require knowledge of the wavelength and 
Ao at several upstream locations so that an estimate 
of the vorticity supplied to the leading edge could be 
made. 

Figure 8 shows the spanwise variation of Frossling 
number at 0 = 0 for the same cases shown in Figure 
7. Notice that locally, the Frossling number drops be- 
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Figure 8: Frossling Number vs. Spanwise Distance, 
Ao — 0.04, Rcr = 6500 

low the two-dimensional result. This reduction corre- 
sponds to the region directly in front of the lower speed 
flow and beneath the vortex pair. The drop in Fr is 
more than offset by the increase in Fr over most of the 
wavelength, causing the spanwise-averaged Frossling 
number to be higher than the two-dimensional result, 
as was shown in Figure 7. The implication of this ob- 
servation is that the spanwise-averaged temperature 
will be lower, but locally the peak temperature will be 
greater than the two-dimensional value. 

To demonstrate the grid dependence of the present 
calculation, Figures 9 and 10 show the Frossling num- 
ber at the leading edge as a function of spanwise dis- 
tance. Figure 9 shows the results for A = 0.47? and 
Aq = 0.04 for two grids. The difference between the 
two grids is that one has geometric stretching in the 
spanwise direction, which clusters points in the re- 
gion where the vortex forms, while the other grid has 
uniform spacing in the spanwise direction. The mini- 
mum spanwise spacing in the stretched grid is roughly 
a factor of 3 smaller than in the grid with uniform 
spacing in the spanwise direction. Notice that the 
Frossling number is only slightly affected by the change 
in grid. Figure 10 shows the results for A = 0.27? and 
Aq = 0.066. The fine grid, which has twice as many 
grid points in each direction, falls within five percent 
of the normal grid calculation. It is expected that at 
lower values of i4o, the agreement would be even bet- 
ter. 

As was mentioned previously, the spanwise- 
averaged Frossling number was higher than the two 
dimensional result for every case considered. Figure 
11 clearly demonstrates that a small spanwise varia- 
tion in velocity is capable of causing a large increase 



Figure 9: Comparison of Grid Dependence 


in the heat transfer coefficient. It figure also demon- 
strates the minor effect that grid refinement has on the 
results. 

VI Conclusions 

From the present calculations it is hypothesised 
that a minimum level of vorticity must be supplied to 
the leading edge for a vortex to form. Using a simple 
parallel flow analysis to predict the decay of vorticity, 
it is further proposed that the appearance of a station- 
ary vortex structure can be predicted from upstream 
conditions. If a disturbance, dominated by a single 
wavelength, is produced upstream of a leading edge, 
then the maximum distance between the disturbance 
and leading edge that will allow the formation of a 
stationary vortex structure will vary linearly with the 
Reynolds number based on the dominant wavelength. 
This statement is shown to agree qualitatively with 
experimental observations [4, 3] as well as with the 
present numerical solution of the Navier-Stokes equa- 
tions. 

The parallel flow analysis also indicates that an 
optimum wavelength exists for a given distance to the 
upstream disturbance and a given Reynolds number 
based on that distance. This optimum wavelength is 
found to vary as the inverse of \/Re^, This result 
may explain the appearance of spanwise variations in 
experiments, with wavelengths that vary inversely as 
the square root of the Reynolds number, as opposed 
to being explained by an instability in the stagnation 
flow as has been proposed in the past. If upstream 
of the leading edge, a spanwise variation exists con- 
sisting of many different wavelengths, then the wave- 
length supplying the most vorticity to the leading edge 
would vary as the inverse of from the above 
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Spanwise Distance, z/X 


Figure 10: Comparison of Grid Dependence 

argument. If Reji varies while holding the ratio Ax/7i 
fixed, this would imply a variation with which 

is what has been seen experimentally. It should be 
noted that even if the decay is not occurring exactly 
as the present analysis would suggest, the trends may 
still remain unchanged. For instance, Townsend [32] 
has indicated that the amplitude of a periodic distur- 
bance in mean velocity for turbulent flow should decay 
as in equation 7, except that the coefficient in the ex- 
ponential would be increased by the ratio of the turbu- 
lent to kinematic viscosity. This technique would not 
contradict any of the above mentioned trends. 

For all cases considered in this study, it was 
found that introducing a spanwise variation into the 
freestream always caused an increase in the spanwise- 
averaged heat transfer coefficient. The percentage in- 
crease in the heat transfer coefficient was substantially 
larger than the freestream disturbance expressed as a 
percentage of the freestream velocity. For example, a 
0.04 disturbance with wavelength 0.4i? located 9 radii 
upstream of the leading edge resulted in an increase 
in the heat transfer coefficient of 18% above the two 
dimensional result. This result demonstrates that a 
small spanwise variation upstream of a body can be 
an important factor in determining the heat transfer 
to that body. 
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